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Lagrangian perturbation theory for cosmological fluid describes structure formation in the quasi- 
nonlinear stage well. In a previous paper, we presented a third-order perturbative equation for 
Lagrangian perturbation with pressure. There we considered only the longitudinal modes for the 
first-order perturbation. In this paper, we generalize the perturbation, i.e., we consider both the 
longitudinal and the transverse modes for the first-order perturbation. Then we derive third-order 
perturbative equations and solutions. 
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I. INTRODUCTION 



The structure formation scenario based on gravitational instability has been studied for a long time. The Lagrangian 
perturbative method for the cosmological fluid describes the nonlinear evolution of density fluctuation rather well. 
Zcl'dovich 1] proposed a linear Lagrangian approximation for dust fluid. This approximation is called the Zel'dovich 

IS EE USB III- ZA describes the evolution of density fluctuation better than the Eulerian 
121 . Once ZA established, the second- and the third-order perturbative solutions for dust fluid 
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were derived [TMlMlll 13 111 E S • 

Recently the effect of the pressure in the cosmological fluid has been considered. Buchert and Dommguez |2l| showed 
that when the velocity dispersion is regarded as small and isotropic it produces effective "pressure" or viscosity terms. 
Furthermore, they posited the relation between mass density p and pressure P, i.e., an "equation of state." Recently, 
Buchert and Dommguez (2^ discussed the origin of "adhesion" in the structure formation scenario. They considered 
two cases, velocity dispersion and coarse-graining. They named these models the Euler- Jeans-Newton (EJN) model 
and the Small-Size Expansion (SSE) model, respectively. 

In this paper, we treat only the EJN model. Following our previous paper |23| . we call this model the "pressure 
model." Adler and Buchert |2 1| have formulated the Lagrangian perturbation theory for a barotropic fluid. Morita 
and Tatekawa |25| and Tatekawa et al. |20l i solved the Lagrangian perturbation equations for a polytropic fluid up to 
the second order. In the aforementioned paper |23j . we solved third-order perturbation equations in the simple case. 

In this paper, we notice the effect of the primordial vorticity. In dust models, Buchert [lfj investigated the behavior 
of the vorticity with the first-order Lagrangina perturbative equations. He showed that the primordial vorticity is 
amplified in proportion to the enhancement of the density fluctuations by deriving the relation between the vorticity 
and the density fluctuations. Barrow and Saich [13j showed the effect of the vorticity in gravitational collapse in an 
expanding Universe. Sasaki and Kasai [2(]] derived the solutions for the perturbations of both the longitudinal and the 
transverse modes up to third-order in E-dS Universe model. Then they showed how the vorticity affects the evolution 
of the density fluctuation. 

In the pressure model, although we showed the perturbative equations for a polytropic fluid up to the second order 
[25). we ignored the effect of the first-order perturbation in the transverse mode. As showed in the past paper |2C| . 
the effect of the vorticity becomes important in nonlinear evolution of the density fluctuation in the dust model. 
In the pressure model, the effect of the vorticity would becomes important. Because the evolution of the first-order 
perturbation in the transverse mode is independent of the scale of the fluctuation, the effect especially appears in small 
scale where the longitudinal mode shows decaying oscillation. Although we presented the third-order perturbative 
equation for the pressure model |23j |. we ignored the effect of the primordial vorticity. The transverse modes do 
not have a growing solution in the first order, we could neglect the transverse modes for the linear approximation. 
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However, the combination of longitudinal and transverse modes do grow in third-order perturbation. Therefore when 
we observe the transverse mode in the initial condition, we need to estimate the effect of the transverse mode. Then 
we derive the perturbative solution for the simplest case, i.e., where the background is given by the Einstein-de Sitter 
(E-dS) Universe model and the polytropic index 7 = 4/3. 

This paper is organized as follows. In Sec.^ we present a Lagrangian description of basic equations for cosmological 
fluids. In Sec. II I II we show perturbative equations and derive perturbative solutions for the pressure model up to 
the third order. In Sec. IIII Al we show the perturbative equations and solutions for generic polytropic fluid in the 
E-dS Universe model. In Sec. IIII Bl and Sec. IIII CI we show the perturbative equations for second- and third-order 
perturbation, respectively. In Sec. II VI we derive perturbative solutions. In general, it is extremely hard to solve higher- 
order perturbative equations. Therefore we show generic solutions for first-order perturbation only. In Sec. II V Al and 
Sec. lIVBl we derive the perturbative solutions for second- and third-order perturbation, respectively. Here we derive 
the solutions only for the simplest case (E-dS Universe model, 7 = 4/3). In Sec.[V] we summarize our conclusions. 



II. BASIC EQUATIONS 

We briefly introduce the Lagran gian description for cosmological fluid. The basic equations for cosmological fluid 
were showed in the previous paper [23l |. Here we show only the evolution equation for Lagrangian perturbation. 

In Newtonian cosmology, to introduce cosmic expansion, we adopt the coordinate transformation from physical 
coordinates to comoving coordinates. 

T 

x = — 77T, a(t) : scale factor, (1) 
a(t) 

where r and x are physical coordinates and comoving coordinates, respectively. The scale factor satisfies the Fried- 
mann equations. 

In the Lagrangian perturbation theory, instead of density fluctuation the displacement from homogeneous distribu- 
tion is regarded as a perturbative quantity. 

r = ax = a(q + s(q, t)) , (2) 

where x and q are the comoving Eulerian coordinates and the Lagrangian coordinates, respectively, s is the displace- 
ment vector. In the Lagrangian model, the displacement vector is regarded as a perturbative quantity. From Eq. J5J, 
we can solve the continuous equation exactly. Then the density fluctuation is given in the formally exact form. 

'- 1 - rl - J - d -(^)- (3 » 

J means the Jacobian of the coordinate transformation from Eulerian x to Lagrangian q. Therefore, when we derive 
the solutions for s, we can know the evolution of the density fluctuation. 
The peculiar velocity is given by 

v = as . (4) 

Then we introduce the Lagrangian time derivative: 

d d 1 _ 

at at a 

Taking divergence and curl of Euler's equation, we obtain the evolution equations for the Lagrangian displacement: 

V,. [s + 2-s~^—J-^V x j\ = -AnGptiJ- 1 -!), (6) 



a 



V I x(H2-s] = 0, (7) 



where (') means the Lagrangian time derivative (Eq. 10). To solve the perturbative equations, we decompose the 
Lagrangian perturbation to the longitudinal and transverse modes: 

s = V5 + s T , (8) 
V-s T = 0, (9) 
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where V means the Lagrangian spacial derivative. 
Here we expand the Jacobian: 



J 1 ~l~ ^i.i ~t~ (^i^iSjJ Sj^j *>j.i) ~t~ det (s^j) 



= i + v 2 5 + ~ {(V 2 S) 2 - S tij S di - s^sli - 23,^} 

+ det(S, ij +sl j ) . (10) 
Because Eqs. © and (J7J) include the Eulerian spacial derivative, we change to the Lagrangian spacial derivative. 
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III. THE LAGRANGIAN PERTURB ATIVE EQUATIONS 
A. The first-order perturbations 

From Eqs. (JHJ and 0, we obtain the first-order perturbative equations: 

V 2 ( S« + 2-SW - tatGpbSW - ^^V 2 S^) = , (11) 
\ a a z J 

a 



V x [m t ^+2-s t ^) = 0. (12) 



The hrst-order solutions for the longitudinal mode depend on spacial scale. Therefore the solutions are described with a 
Lagrangian wavenumber. In this paper, we discuss only perturbative solutions in the E-dS Universe model |25il2fill2^ |. 

SW(K,t) = C + (K)D+{K,t) + C-{K)D-(K,t), (13) 



£-1/6^(^^-7+4/3) fo r 7 ^ 
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D ± {K 1 t) - I t -l/6± N /25/36-S|JCP fo r 7= | j ^ 

3y^7(67rG)( 1 ^)/ 2 _ 4 /3 
A= jj-^j , B =r (6nG) , 

where J means Bessel function. The coefficient v is given by v = 5/(8 — 67). C ± {K) is given by the initial condition. 
For the transverse mode, the solutions are same as for dust model: 

S T(1) = C T+ +C T- t -l/3_ (15) 

The transverse mode does not have a growing solution. Therefore, when we consider only large-scale fluctuation, the 
longitudinal mode dominates during evolution. 

According to a comparison of the first-order perturbation with a numerical simulation |28j . the first-order pertur- 
bation seems rather good until the quasi-nonlinear regime (5 ~ 1). After that, we need to consider higher-order 
perturbation. 

B. The second-order perturbative equations 



From Eqs. ||BJ and Q, we obtain the second-order perturbative equations |25L l2T| . For the longitudinal mode, the 
equation becomes 

V 2 [ + 2^ 2 > - 4nG Pb S^ - ^Llv 2 S( 2 M = Q L ^ , (16) 
\ a ar I 
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Q L ^ = 2nGp b 



7-I 



s$}s® - (v 2 5 (1) ) ; 

{ 7 V 2 5j ) V 2 5y ) + (7 - 1)V 2 V 2 S«V 2 S« + + 2S^W 2 S^} 



7-1 



{sgv 2 ^ 



27rr^s T(1) s T(1) K7P b 1 
-zirupbs l j a ii ^ — 



T{1) T(l)"l 



For the transverse mode, after some arrangement, we can describe as follows: 



V 2 



T(2) 



(17) 



(18) 



Qf (2) - -^{^^g + ^fV 2 ^-^^ 1 ^^-^^} 

{sgLv's$ + s^sf - V 2 ^V 2 sg - ^V 2 S« } 

+ (-JX^ + ^Vfl? - W^W^S - sj^S%) . (19) 

Here we can note second-order transverse mode solutions. In the pressure model, even if we consider only the 
longitudinal mode for the first order, the second-order perturbation for the transverse mode appears. In dust model, 
it does not appear. Therefore, when we derive the third-order perturbative solutions, we must consider the second- 
order transverse mode. 



C. The third-order perturbative equations 



In the previous paper, we show the third-order perturbative equation in which only the longitudinal modes for the 
first-order perturbation considered [2]|. The third-order perturbative equation becomes very complicated. Here we 
introduce scalar quantities generated by the Lagrangian perturbations. 



/*!(*) EE V 2 *, 

H 2 (A,B) ee ~((V-A)(V-B)-Ai tj B jti ) 

p, 2 (A) ee p 2 (A 7 A), 

p 3 (A) ee det(A (J ) , 



(20) 

(21) 
(22) 
(23) 



where A and B are vector quantities. ^ is a scalar quantity. Using these quantities, the Jacobian (Eq. (|10f> ) is written 
as 



J = pi(S) + p 2 {S,t + sf) + p 3 {S,t + sj) . 



(24) 



First, we show the longitudinal mode equation. As in the second-order perturbative equation, we separate the 
terms of the third-order perturbation from the others. Then the terms consisting of the first- and the second-order 
perturbations are collected to the source term Q L ^\ 



.7-1 



(25) 



We consider the source term. Using Eqs. (|2U[I - H23[I . the terms are written as follows: 
Q m = 4nGp b \p 1 (S^f+p 2 (S? +S ^ 1 \s^+ S ^ 
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+ (^ + «D 1, )(^ , +^;j +J ™ + 2^') 
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T(2), 



?(!)> 



-{(5i i 2+-a i) )(/*i(5 (a) )+/* 9 (si l) +«r (l) )), 
-{ 7M1 (^))( M1 (^))+ M2 (^+.sr 

{7 (mi(^ (2) ) + + *r (1) )) V 2 ^ - ^f±W (1) ) 2 MS (1) h) 

+7{Mi(S (1) )(S < J 1 2+-a i) )/'i(S (1) ),*} J 

+(4 1 2+-3 1 M^V+-3 1) )A*i(^)! I } ifc 

+7(^+ S 3 1) ){Mi(^ 1) W(^ 1) ).} fc 



(26) 



The transverse mode also seems complicated. However, if we neglect the first-order transverse mode, the evolution 
equation is described as 



(C (3) + 2^f (3) )=Qf (3 \ 



(27) 



T(2) 



,k 



-Us 



, r(i) 



o(2) , 9 a «j(2) 

a ,ij 



^) +2 ^n 2) 



t,3 



a 



+ < S 



^G Pb sf k 



'2°) 



+ 



7-1 
«7P b 



V/2 c(l) c(l) 



W 2 c(1)t7 2 q« , V7 2 C« cC 1 ) qW 
v ° : ife v °j D ,jk v ° ,ikl° ,jk° ,jl 



c(1)v72q(1) c(1)c.(1) V7 2 q( 1 ) 
' D ,ijk D , 3 l V D ,kl D JfcZ V 



+47rGp 6 



T(l)o(l)c(l) 
S Mj D ,fcZ D ,jk 



, T(l) / 2 „(!)„(!) V 2c(l)a(l) , o(l)c(l) rfloW 1 cWcW , oWrf)) 
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2 T(l) T(l) 



r(i) r(i)\ <,(i) 



T(l) T(l) 
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(T(l) T(l) T(l) T(l)\ 

+(v 2 s r (i) ^ ) +s s i) ^ ) ) v25 



(i) 



r(i) r(i) v2<? (i) 

S l,j S k,l V D ,r/fc 



(28) 



IV. LAGRANGIAN PERTURB ATIVE SOLUTIONS FOR HIGHER-ORDER PERTURBATION 
A. The second-order perturbative solutions: The simplest case 

The second-order solutions are formally written as follows: 

= 



^ J &t'G{K,t,t')Q L( ?\K,t) : 
dt'G T (t, t')Q T ^{K,t), 



\K 
1 

G T (t,t') = 3(i' - t-VV 4 / 3 ) . 
G L (K, t, t 1 ) depends on the "equation of state." If 7 7^ 4/3 and v — 5/(8 — 67) is not an integer, we have 



G L {K,t,t') 



tx 4 

: -7+ - 

2 sin vn \ 3 



-Ju (A\K\r 



t -i/6 t n/e 
" 1 ' ''\j_ v (A\K\t'-~< +i / 3 



J- v (A\K\t-~<+^) J v (A\K\t>- 



7+4/3 



and if 7 = 4/3, 



G L (K,t,t') = --*t- l '*ir'*(rWit*V-&v 1 >-KKj} , 



UK) ee yj 2 ±-B\K\\ 



(29) 

(30) 
(31) 



(32) 

(33) 
(34) 



The second-order perturbative solutions for the case of 7 = 4/3 in the E-dS universe model has already been derived 
by Morita and Tatekawa (25J. The solutions are described by 



dK' (3i(K, K')Ei(K, K', t) + E^(K, K')E^(K, K', t) 



S T ( 2 ) 



(K) = J dK' (Ej^K^^EfiK^K^^+ElyK^K'jE^iK^K'y) 



(35) 
(36) 



{{K 1 ■ (K - K')) 2 - \K'\ 2 \K - K'\ 2 } 



-l\2 



(2tt) 3 \K\ 2 

+b\\\K'\ 2 \K - K'\ 2 {K' ■ (K - K')) + \\K'\ A \K - K' 

+ {{K 1 ■ (K - K')f + 2\K - K'\ 2 {K' ■ (K - K')) 2 ) }] 



iB 1 

WW 



K'\ 2 {K' ■ (K - K 1 )) + 2(K' ■ (K - K 1 )) 



l\\2 



(37) 
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+ \K'\ 2 \K - K'\ 2 + 2\K'\ 2 (K' ■ (K - K'))\ K\ , 



1 1 

i B 



~(K' ■ (K - K')) -B{\K- K'\ 2 + (K' ■ (K - K'))} 
K - K'\ 2 {K' ■ (K - K 1 )) 



(38) 

K[{K 3 -K'), (39) 



(2tt) 3 \K\ 2 ' 

[(K> ■ (K K'))K[ + \K- K'\ 2 K[ \K'\ 2 {K, - K[) - (K> ■ (K *"))(#, %)] , (40) 
{{{K'-(K-K')) + \K-K'\ 2 )K' t 



(2tt) 3 



- (|Kf + (K> ■ (K K'))) {Ki - Ki)} (K 3 K' s ) 
-B{((K' ■ (K - K'))\K - K'\ 2 + \K- K'\ 4 ) K[ 

- (\K'\ 2 \K K'\ 2 -\K- K'\ 2 {K' ■ (K K'))) {K t K<)} (K, - ty] 



(41) 



Ef(K,K',t) cx 



E^(K,K',t) cx 



E^AK,K',t) cx 



2b(K) 



E E 

(e=±) (®=±) 

E E 



C®(K')C®(K - if')i-i/3©6(-Rr')®6(K--Jc') 
((®b(K') i 6(K - JsT') - i) 2 - 6(K) 2 ) 
C r ®(K" / )C" 8 (if - K')t- 1 /3©i/6®6(J<:') 



C T+ {K')C T+ {K - K') ( C^+^QC 71 -^-^) _ 1/3 



26(AT) 



3- - b(iiT) 2 



+ 



36 

cF-iic^+^K - K') t _ 1/3 C T -(K')C T -(K-K') t _ 2/3 



h-HK) 2 



J_ 

Mi 



El{K,K',t) cx 3^ 

(®=±) (®=±) 



- &(*0 2 3 " HK) 2 

C®(K')C®(K - K')t- 1 / 3 ' sb( - K ")® b( ' K - K ") 
((Bb(K') b(K - K') - \){®b{K') ® 6(Jf - K')) 



E^{K,K',t) cx £ 



(e=±) 



where Yl(®=±) means 



C T ®{K')C+{K-K>) 1/6mK) 
KK) 2 -± 

C T ®(K')C~(K - K') 

h (b(K - K>) ® (-1)) (b(K K>) © (-1)) ' 



^ (a® © /?) = (a+ + /?) + (a" - 0) . 
(e=±) 



-i/2eb(-ff--fO 



(42) 
(43) 



(44) 
(45) 



(46) 



(47) 



E% and -B^ appears only in the pressure model. If we ignore the effect of the pressure, only E±, E^, and E% appear. 
From these results, it is apparent that the pressure does affect the evolution of the higher-order perturbation. 



B. The third-order perturbative solutions: The simplest case 



As in the first- and the second-order solutions, the third-order solutions are described with the Lagrangian wavenum- 
ber. Following the method used in the second-order solutions, the third-order solution is given by this integration: 



S {3) = -j^ £ dt'G(K,t,t')Q L &(K,t) , 



(48) 
(49) 



In the previous paper [23|], we showed the third-order perturbative solutions for the simplest case, the case of 7 = 4/3 
in the E-dS Universe model. In this case, the contributions of the gravitational terms and the pressure terms become 
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identical: 

7-1 

AnGp b , — cx a- 3 . (50) 

In addition to these conditions, we considered only the longitudinal modes in the first-order perturbation. Here we 
consider both the longitudinal and the transverse modes in the first-order perturbation. As we show in Sec. lIIICl the 
third-order perturbative equations become more complicated. To avoid complexity, we estimate only time evolution 
in the third-order perturbative solutions F(t): 

F(t) cx t a . (51) 

The third-order perturbation is affected from the first- and the second-order perturbation via source terms Q L and 
Q T : 

gL(3) = Qm ^1)^ S {2) ; 5 <1) „ s T<a) )S ni)* 5 (2) jS T(l)^ s T(2) ) 

S« * S« * S (1 \S^ * SW * s T ^,S^ * s T « * s T « s T « * s T ^ * s T ^) , (52) 

Q T(3) = QT(3) ^(l)^(2) )5 (l)^ s T(2) )S r(l)^( 2 ) )S r(l)^ s T(2) ) 

* <?« * S« S« * S« * s T « S« * S T « * S T «) , (53) 

where * means convolution. 

In Sec. IIII Al and IIV Al we showed the time dependence of the perturbative solutions. Here we arrange and list 
again: 

D ± {K) oc r 1 / 6 ^*), (54) 

D T± CX ^-l/6±l/6 ; (55) 

Ei(K, K') cx r V3±6(*')±Kx-if') , (56) 

E%(K,K') cx r V3±i/6±6(*-JO ) (57) 

E£(K,K') cx r V3±i/6±i/6 ? (58) 

Ef(K,K') cx r V3±6(if')±Kx-if') , (59) 

£ 2 T 0fi:,ii:') cx r i/3±V6±6(*-K') . (60) 

From Eq. I|34|l . if we take a limit of weak pressure (P — > 0), converges to 

W-g- (61) 

Here we notice the relation between i?f and , i.e., time evolution in both the longitudinal and the transverse 
modes. Both of these terms are generated from * S^- 1 ' in Q L ^ and Q T ( 2 \ Therefore if the convolution in the 
source terms in both modes are identical, the evolutions will also be identical. 

Tables HI and UTI show the time dependence of the third-order perturbative solutions. When the dust model |20( and 
the pressure model are compared, we notice that several terms in the third-order perturbative solutions do not appear 
in the dust model. For example, the contribution of the term DEf , which is proportional to DE±, i.e., one of most 
dominant term, appears only in the pressure model. 



V. SUMMARY AND CONCLUDING REMARKS 



In this paper, we showed the third-order Lagrangian perturbative equations for the cosmological fluid with pressure. 
Then we derived third-order perturbative solutions for a simple case. 

In the pressure model, we must notice one other point. The longitudinal mode in the first-order perturbation 
depends on the scale (or the wavenumber) of the fluctuation. On the other hand, the transverse mode in the first- 
order perturbation is independent of the scale of the fluctuation. In the dust model, because the transverse mode 
does not have growing modes, even if the primordial vorticity exists, we can simply ignore this mode. However in 
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TABLE I: The time dependence of the third-order perturbative solutions in the longitudinal modes (the case of 7 = 4/3 in 
the E-dS universe model) . The contribution means the time component in the source terms. The indices of the time evolution 
a (Eq. 15110 depend on the convolution in the source terms Q L ^ 3 \ Here we define b 1 = b(K"),b2 = b(K' — K"),bz = 
b(K — K'), 64 = b(K — K"). "P" denotes that the contribution does not appear in the dust model. 



The convolution 




The contribution 


The indices a 




S w (K")*S i2 \K - 


K") 




-| ±61 ±b 2 ±b 3 








DE2 


-|±I± 6l ±6 4 


P 






DE^ 


-|±|±|± 6l 




S w (K")*s T ^(K - 


- K") 


deT 

DEI 


-|±6i ±b 2 ±h 
-|±i±6i±6 4 


P 


T(l) 1 rs-tt\ a(2) { jv- 

s K } (K ) * S y } {K - 


- K ) 


D T Ei 


-| ± | ±6 2 ±63 








D T E2 


-!±!±!±&4 


P 






D T Et 


-i ± i ± i ± i 




s T ^(K")*s T ^(K 


— K") 


D T El 
D T El 


-i±i±6 2 ±63 
-|±I±I± 64 


P 


S m (K")*S w (K' - 


-K")*S W {K -K') 


DDD 


-i±6i±6 2 ±b 3 




s T ^(K")*S w (K' 


-K")*S W (K -K') 


D T DD 


-I ±| ±6 2 ±63 






-K")*S W (K -K') 


D T D T D 




P 


s T ^(K")*s T ^(K' 


-K")*s T{1 \K-K') 


D T D T D T 


2 - 1 - 6 6 - 1 - 6 





TABLE II: The time dependence of the third-order perturbative solutions in the transverse modes (the case of 7 = 4/3 in 
the E-dS universe model). The contribution means the time component in the source terms. The indices of the time evolution 
a (Eq. 1511 ) depend on the convolution in the source terms Q L ^\ Here we define b\ = b(K"),b2 = b(K' — K"),bs = 
b(K — K'), 64 = b(K — K"). "P" denotes that the contribution does not appear in the dust model. 
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The indices a 
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DEt 
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1 
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± 61 ± 62 ± 63 
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P 
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D T Et 
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the pressure model, we cannot ignore it. The growing factor in the longitudinal mode depends on the scale of the 
fluctuation; when we consider small-scale fluctuation, the transverse mode dominates. 

The transverse modes do not appear in spacial one-dimensional model. In two- or three-dimensional model, the 
computation of the perturbation becomes huge. Especially, in the third-order perturbation, enormous mode-coupling 
computation is required. Even if in one-dimensional model |23| . the computation of the third-order approximation 
required long time. When we compute the third-order approximation in generic case, we will apply Monte Carlo 
integration .29] for the computation of the mode-coupling. 

Here we estimate importance of the transverse mode in the pressure model. For simplicity, we consider only the E-dS 
Universe model. In past papers |25ll26|. we showed that when the polytropic index 7 is less than 4/3, all longitudinal 
modes oscillate and decay. On the other hand, when the polytropic index 7 is more than 4/3, all longitudinal modes 
grow as those in the dust model do. The case of 7 = 4/3 is critical. We consider the index of the time component in 
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the first-order perturbation (Eq. (|14[l h If the Lagrangian wavenumber K satisfies the condition 



V oB 5 

the longitudinal modes dominate during evolution. Therefore when we consider small-scale (\K\ > (2^/6/5)Kj) 
fluctuation with primordial vorticity, we must consider the evolution of the transverse modes. 

Recently several dark matter models have been proposed [30| . In several model, special interaction affects avoidance 
of matter concentration. In general, these interaction is introduced by scattering cross-section. Instead of the cross- 
section, if the interaction in some kind of dark matter can be described by the effective pressure, we can examine 
the behavior of the density fluctuation in a quasi-nonlinear stage. Furthermore, when we compare the observations 
and the structure that is formed by using the pressure model, we can delimit the nature of the dark matter. In 
future, when we consider the evolution of the density fluctuation or the peculiar velocity in even high-density regions, 
although these seem highly complicated, the third-order perturbative solutions may become useful. 
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